Reproduced  From 
Best  Available  Copy 


College  of  Engineering 

Virginia  Polytechnic  Institute  &  State  University 
Blacksburg,  Virginia  24061 


VPI-E-74-10 


A  FINITE  DEFORMATION  ANALYSIS  OF  THE 
NEAR  FIELD  SURROUNDING  THE  TIP  OF 
CRACK-LIKE  ELLIPSES 

J.  J.  McGowan*  and  C.  W.  Smith** 


Department  of  Engineering  Science  &  Mechanics 


Reporting  Period: 


May  1974 


i- TO-  /-r?' 

R » <■  Uii, .luiC}}}  4 r 01  z’  .  T: ,  T  - 


^  "  Tu  ;  U 


1973-74  Academic  Year 


affdl/for 

V)/PAFB,  OHIO  45433 


Prepared  For: 


National  Science  Foundation 
Washington,  D.  C. 

Grant  No.  GK  39922 


Approved  for  Public  Release;  distribution  unlimited. 


20000111  065 

*Assistant  Professor 

**Professor  - 

Department  of  Engineering  Science  &  Mechanics 
Virginia  Polytechnic  Institute  &  State  University 
Blacksburg,  Virginia  24061 


1.  Report  No.  2.  Government  Accession  No. 

VPI-E-74-10 

3.  Recipient's  Catalog  No. 

4.  Title 

5.  Report  Date 

A  FINITE  DEFORMATION  ANALYSIS  OF  THE  NEAR  FIELD 
SURROUNDING  THE  TIP  OF  CRACK-LIKE  ELLIPSES 

May  1974 

8.  Performing  Organization 

7.  Authors 

Report  Number 

J.  0.  McGowan  and  C.  W.  Smith 

VPI-E-74-\0 

9.  Performing  Organization  Name  &  Address 

10.  Work  Unit  No. 

Department  of  Engineering  Science  &  Mechanics 

Virginia  Polytechnic  Institute  &  State  University 

VPI -32371 4-1 

Blacksburg,  Virginia  24061 

11 .  Contract  or  Grant  No. 

12.  Sponsoring  Agency  Name  and  Address 

GK  39922 

National  Science  Foundation 

13.  Type  of  Report  & 

Washington,  D.  C. 

16.  Abstract 

Period  Covered 

Academic  Year 

1973-74 

A  finite  deformation  analysis  of  the  region  surrounding  the  tip  of  small 
crack-like  ellipses  in  an  infinite  plate  under  all  around  tension  is  presented. 
The  study  is  carried  out  in  the  deformed  geometry  and  includes  the  effects  of 
finite  strains  and  rotations.  A  stress  function  is  first  introduced  to  the 
complete  compatibility  equations  through  linear  constitutive  relations;  the 
resulting  governing  equation  is  solved  through  finite  differences.  The  range 
of  root  radii  investigated  varies  from  one  to  nine  times  that  of  a  deformed 
crack.  ^Normal  stresses  and  strains,  maximum  in-plane  shear  stress,  displacements 
and  stress  intensity  factors  are  presented.  The  results  of  the  analysis  are 
compared  to  the  linear  analysis  of  Inglis.  The  effects  of  finite  strains  and 
rotations  are  shown  to  be  large  but  are  concentrated  within  a  few  root  radii 
of  the  tip.  lit  has  been  found  that  two  types  of  behavior  exist  at  the  tip  of 
notches,  depending  on  the  size  of  the  root  radius:  small  root  radii  produce 
maximum  stresses  away  from  the  tip,  whereas  larger  root  radii  produce  maximum 
stresses  at  the  tip  itself,  'i 


17.  Key  Words 

Stress  Fields,  Crack  Tips, 
Finite  Deformations, 

Fracture  Mechanics 

19.  Security  Classif.  (report) 
Unclassified 


18.  Distribution  Statement 

Approved  to  Public 
Release  Distribution 
Unlimi  ted 

20.  Security  Classif.  (page) 
Unclassified 


NOMENCLATURE 


Coefficients  as  defined  in  (2.7) 

Base  vectors  in  the  deformed  system 
Dimensional  1/2  crack  length  as  shown  in 
Figure  2.2,  mm 

2 

Young’s  modulus,  N/mn 
Permutation  tensor  defined  by  (2.6) 
Coefficients  as  defined  in  (2.5) 

Coefficients  as  defined  in  (3.1) 

Covariant  and  contravariant,  respectively, 
deformed  metric  tensors  nondimensionalized  by 

9ij  "  9ij/^  ’ 

Covariant  and  contravariant,  respectively, 
undeformed  metric  tensors  nondimensionalized 
by  hij  =  h^j/F^  ;  h^^  = 

Stress  intensity  factor,  N/mm^^^ 

Theoretical  stress  intensity  factor,  N/mm^^^ 
Apparent  stress  intensity  factor,  N/mm^^^ 
Position  vector  in  the  deformed  and  undeformed 
systems,  respectively,  mm 
Covariant  and  contravariant  displacements, 
respectively,  nondimensionalized  by  w^-  =  w^./F 
=  w’/F 


iv 


X,  y 

Local  cartesian  coordinates  nondimensionalized 

by  X  =  x/c,  y  =  y/c 

Coefficients  as  defined  in  (3.2) 

A  ,  A 1  3  A2  5 

Coefficients  as  defined  in  (3.1) 

e 

Perturbation  parameter 

Covariant  strain  tensor 

^ijkl 

Strain  functions  as  defined  by  (2.1) 

-> 

c 

Displacement  vector 

n  >  C 

Local  orthogonal  curvilinear  coordinates  as 

shov/n  in  Figure  2.2 

0 

Sum  of  physical  normal  stresses 

V 

1 

Poissons  ratio 

^0 

Value  of  ?  at  boundary  of  ellipse 

p 

Ellipse  root  radius  nondimensionalized  by 

p  =  p/c 

0 

Dimensional  remote  normal  stress  3  psi 

aij  3 

Covariant  and  contravariant,  respectively 3 

stress  tensor  nondimensionalized  by 

"^max 

Maximum  in-plane  shear  stress  nondimension¬ 
alized  by  “  '’^max/'^ 

$  5  (p 

Outer  and  inner3  respectively,  stress  function 

nondimensionalized  by  $  =  (j)  = 

1A| 

Determinant  of  A-jj 

V 


Superscri pts 
Dimensional  quantity 
Physical  quantity 
Inner  variable 
Outer  variable 
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I.  INTRODUCTION  AND  REVIEW  OF  THE  LITERATURE 


Although  the  use  of  photoelasticity  in  obtaining  stress  intensity 
factors  was  noted  by  Irwin  [1]  in  the  early  1950's,  broad  implementa¬ 
tion  of  the  method  was  slow  to  develop.  The  reason  for  this  is 
believed  due  to  the  problems  involved  in  extracting  accurately  the 
stress  intensity  factor  (SIF)  from  photoelastic  data.  Several  dif¬ 
ferent  approaches  have  been  cited  for  achieving  this  purpose  for  two 
dimensional  geometries:  one  developed  by  Marl  off  and  his  associates 
[2]  in  1970;  one  recommended  by  Irwin  (above)  and  applied  by  Smith  and 
Smith  [3-5];  and  those  developed  by  Kobayashi  and  his  associates  [6-9]. 

An  attempt  to  quantify  the  errors  involved  in  the  determination 
of  the  SIF  in  two  dimensions  through  photoelasticity  has  been  pre¬ 
sented  by  Schroedl ,  McGowan  and  Smith  [10].  In  this  study  strong 
emphasis  was  put  on  taking  photoelastic  data  in  a  region  where  Irwin's 
two  degree  of  freedom  system  would  apply  (referred  to  as  a  'singular' 
zone).  Outside  of  this  zone  effects  of  free  surfaces  other  than  the 
crack  surface  may  dominate.  A  Taylor  Series  expansion  approach  (TSCM) 
was  developed  in  [lo]  to  account  for  these  'outer'  effects.  Inside  of 
this  'singular'  zone  crack  tip  blunting,  finite  rotations  and  finite 
strains  dominate  the  desired  crack  behavior.  In  [10]  the  crack  tip 
blunting  is  treated  superficially  and  the  last  two  effects  are  not 
included.  Near  the  tip  of  the  deformed  crack  these  effects  are  very 
significant  and  cannot  be  neglected. 
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In  many  materials  tensorial  nonlinearities,  as  discussed  by 
Ronay  and  Freudenthal  [11]  and  Poynting  [12],  play  a  large  role  when 
strains  are  high.  However,  the  materials  used  in  stress-freezing 
photoelasticity  (Hysol's  CP-4290  and  Photolastic'S  PLM-4B  and  PSM-8, 
for  example)  do  not  exhibit  any  stress-strain  nonlinearities,  even  at 
high  stress  levels,  as  reported  by  Fessler  and  Mansell  [13]  for 
Araldite.  These  tensorial  effects  are,  therefore,  neglected. 

Determination  of  values  of  the  SIF  in  three  dimensional 
problems  is  discussed  by  Smith  [14].  Although  Kassir  and 
Sih  [15]  have  shown  that  locally  along  an  embedded  elliptical  flaw  a 
two  dimensional  representation  is  still  possible,  the  'valid  data'  zone 
may  become  so  constricted  that  the  'inner'  and  'outer'  effects  must  be 
dealt  with.  Smith  et  al  [16-19]  have  developed  the  TSCM  approach  for 
various  geometries,  effectively  handling  the  'outer'  problem.  It 
therefore  becomes  increasingly  important  to  find  a  method  for  accounting 
for  not  only  blunting  but  also  finite  rotations  and  finite  strains  in 
the  tip  region,  the  fracture  aspects  of  which  are  described  in  [20]. 

Several  approaches  have  been  used  for  dealing  with  the  non-linear 
effects  in  the  'inner'  zone.  Perhaps  the  best  known  of  these  is  exempl- 
fied  in  the  solutions  of  Rice  and  Rosengfen  [21]  and  Hutchinson  [22]. 
These  solutions  imposed  small  strain  restrictions  on  an  elasto-plastic 
solution  using  a  detormation  type  plasticity  theory  with  power  law 
constitutive  relations.  From  the  form  of  the  linear  solution  a  class 
of  functions  was  selected  to  represent  the  stresses  and  strains  near 
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the  crack  tip  in  an  asymptotic  sense.  The  solution  was  reduced  to  an 
eigenvalue  problem  which  was  solved  numerically  and  the  subsequent  sin¬ 
gularity  was  evaluated  through  use  of  the  J  integral.  Ref  [21]  was 
carried  out  for  plane  strain  and  Ref  [22]  to  plane  stress. 

A  second  approach  due  to  Knowles  and  Sternberg  [23],  [24]  concerns 
the  same  crack  and  loading  geometry  as  above;  however,  this  latter 
plane  strain  treatment  includes  the  effects  of  finite  geometry  change 
and  finite  strains.  A  non-linear  constitutive  law  is  postulated  for 
large  deformations  using  a  four  parameter  representation  of  the  strain 
energy  density  in  terms  of  the  strain  invariants.  Again  the  asymptotic 
form  of  the  displacements  is  deduced  from  the  linear  problem  and  the 
investigation  reduced  to  an  eigenvalue  problem.  The  solution  contains 
two  arbitrary  constants,  one  of  which  is  obtained  through  a  conservation 
law.  This  approach  was  then  extended  to  clarify  and  improve  the  lowest 
order  asymptotic  solution. 

Prior  to  either  of  the  above  approaches  0.  L.  Bowie  [25]  studied 
the  effect  of  the  small  geometry  change  assumption  on  the  near  field  of 
an  elliptical  hole  using  classical  theory.  He  found  stresses  of  the 
order  of  Young's  Modulus  and  strains  of  order  unity  at  the  notch  tip. 

More  recently,  Mansfield  [26]  employed  a  procedure  involving  an  incre- 
mentalized  classical  theory  of  plane  strain  in  order  to  estimate  stresses 
near  the  tip  of  a  deformed  crack  which  he  took  to  be  elliptical  in  shape. 
His  stress  increments  were  computed  from  the  deformed  geometry  and  also 
indicated  that  stresses  at  the  crack  tip  were  at  the  order  of  Young's 
Modulus.  Although  neither  of  these  two  studies  utilize  large  strain 
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definitions,  they  do  concern  the  very  near  field  behavior  as  will  the 
model  described  in  the  sequel. 

The  present  study  will  focus  on  the  behavior  of  small  crack¬ 
like  ellipses  (the  size  of  these  ellipses  will  be  the  same  as  a  de¬ 
formed  crack)  in  an  infinite  plate  in  biaxial  tension.  This  investi¬ 
gation  seeks  to  give  a  complete  treatment  of  the  results  of  crack-tip 
blunting,  finite  rotations  and  finite  strains  while  still  using 
linearly  elastic  constitutive  relations.  The  problem  will  be  studied 
in  the  deformed  geometry  to  simplify  the  boundary  conditions  and 
overall  approach. 

In  the  first  step  of  the  investigation  a  generalized  'Airy- type' 
stress  function  is  introduced  through  the  linear  constitutive  equations 
to  arrive  at  the  governing  equation.  Next  a  perturbation  analysis  is 
performed  to  characterize  the  form  of  the  equations  in  the  'inner 
nonlinear  region,  and  to  determine  the  matching  conditions  of  the  non¬ 
linear  to  the  linear  solution  at  some  acceptably  far  distance  from  the 
ellipse  tip.  The  nonlinear  compatibility  equation  is  then  solved  in 
the  tip  region  through  use  of  finite  differences  and  collocation. 
Finally  the  complete  strain-displacement  relations  are  integrated  to 
determine  the  initial  geometry.  The  stress  and  strain  distributions 
of  an  initial  crack  will  then  be  ascertained  as  a  limit  to  the  small 


ellipses  analyzed. 


II.  FORMULATION  OF  THE  PROBLEM 


The  characterization  of  small  crack-like  ellipses  in  an  infinite 
plate  under  biaxial  tension,  including  the  effects  of  blunting,  finite 
rotations  and  finite  displacements  will  now  be  presented.  First  the 
governing  equation  applicable  to  the  problem  will  be  derived  and  the 
tensorial  notation  will  be  introduced.  Finally  the  perturbation  con¬ 
siderations  will  be  discussed,  separating  the  inner  from  the  outer 
regions  (i.e.,  nonlinear  from  linear);  and  the  resulting  solution 
procedure  will  be  outlined,  setting  the  stage  for  the  actual  solution 
method. 


A.  Notation  and  Formulae 


The  generalized  curvilinear  orthogonal  coordinate  system  used 
in  this  study  is  shown  in  Figure  2.1,  and  the  compatibility  equation 
as  given  by  Sokol nikoff  [27]  is 


ijkl 


^jkB 


^ila"  ^jle 


E  -I,  )  ~ 

ika 


(2.1) 


where 

^ijkl  ^jl,ik  ^ik,jl  “  ^jk,il  "  ^il,jk 
"^lik  "  '^ik.j  ^kj,i  "  ^ij,k 

and  h^*^  is  the  contra  variant  metric  tensor  of  the  undeformed  geometry. 
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Note  that  throughout  this  treatment  the  comma  represents  the  covariant 
derivative  with  respect  to  the  deformed  coordinates. 

The  undeformed  covariant  metric  tensor,  h^^^,  as  given  by  [27] 
may  be  v^ritten 


(2.2a) 


or 


f(gn  -  2en) 


“^^12 

(g22  -  2^22) 


(2.2b) 


so  therefore 
hij  = 


r 


I  hi 


(922 

2e 


2622) 


2e 


12 


12 


(Oll  -  2e^^)j 


(2.2c) 


with 


|h|  -  9229-1 1  "  ^9n^22  ~  ^^22^11  ^^22^12  ~  ^^12  * 


The  linear  constitutive  relation  can  be  written  as: 


(1  +  v)  T  _  ye  - 
^jk  "  jk  E  9jk  ’ 


(2.3) 


—  —  i  — ^  -j  — 

where  0  =  0^.  =  g 


i  ,j  =  1 ,2,3. 


For  plane  strain  =  0  if  i  or  j  =  3  so  that  (2.3)  becomes 


"jk 


_  (1  ^  v) 


lUJ-vIbF. 


"jk  ■  E  ®  • 


(2.4) 


-  -IJ  - 

where  0  =  a.  =  a  a_. _. 


i  ,.i  =  1,2. 
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The  only  non-zero,  unique  compatibility  equation  is  for 
ijkl  =  2121.  Noting  this  and  then  introducing  (2.2c)  and  (2.4)  into 
(2.1)  and  collecting  terms  results  in: 
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-  /I  +  v\2i 


(  ^  )  [2{(1  -  v)  022 j  ^11,2 

"1  1 


911  _ 

^  V== -  Or 


,_2)  -  {(1  -  v)  022,1  -  '^“11,1 


-  2<>,2,2lt('  -  '’)  '’,1,2  -  '^22  ‘’22.2  ' 


—  9ll  ._ 

{(1  -  v)  a.|.|  .|  -  v= — '^22,2 
^22 


922  -  . 

■  ,2^  ’ 


h22  = 


1  -  M-  1, 

vj  Qii  -  V —  ^22^  ’ 


1  r77  o  +  ’^)  rn 

iTT  ^  E  l(i 


—  922  _ 

°22  -  1t7'^11^^  ’ 


012)  , 


|hl  =  922  ■  ^  ^22  +  <^22  9n ) 


wl  +  v,2r-  -  /I  \-  -  ^^^22 

+  4(-^)  [on  022  -  ''(1  -  '’)  9n  922  ‘g22  *  911 


-  Oi2  ]  . 


Introducing  now  a  generalized  'Airy'  stress-function  as  sug¬ 
gested  by  Adkins,  Green  and  Shield  [28]: 
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^St  ^Sa 


—  —  -Oil 


^te 


ej  - 

e  $ , 


where 


1 

*111 


0 


i  f  i  =  1 ,  j  =  2 

if  i  =  2,  j  =  1 

if  i  =  j 


Non-dimensional izing  the  following  variables: 

9ll  ^  »  ^22  ^^^22  ’  ^ij  ^ 

=  pa^-j  ,  ¥  =  pc^$  ,  w  =  cw  , 


(2.6) 


where  c  is  some  characteristic  length  and  p  is  some  characteristic 
stress,  and  also  defining  a  physical  derivative: 

Aj  i  =  where  i_  is  not  summed, 

*^gi  i 

and  then  finally  introducing  these  substitutions  plus  (2.6)  into  (2.5) 
leads  to  the  expression: 


o  ~  ^  ^  2 

0  =  {1  +  2eAi  (3>»22  ^’11 )  l-^’ll  ^’22  '  ^’12 

+  ^2  ('J>»22  ^  ^’11^  ^o^^’llll  ^  ^^’1122  ^  ^’2222^ 

+  e{l  +  2e(A3i,ii  +  ^^l^’lll^^  ^2^’122^ 
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+  B3i,in^.T22  ’’’  ^4^’222^’'’  ®5®’112^’'’  ^6^’112^’222^ 

+  e{l  +  2e(A3$j22  '*’  ^^1^*222  ^^^*112 

+  B3i,n2^’222  ^  ^4^’lll^'^  ^5^’122^’^  ®6^’111^’122^ 


+  e^2i,i2  ^22^*112  '*’  ^2^’in^’222 


+  C3(i,^22^’222  ^’111^*112^^ 


(2.7) 


where 


=  f  (1  +  v),  =  -(1  -  2v),  A^  =  v(l  -  v) 


A3  =  -(1  -  v),  A^  =  V,  Bq  =  (1  -  v),  B-j  =  -v(l  -  v) 

B2  =  (1  -  v){2  -  v),  B3  =  1  -  4v  +  2v^,  B^  =  (1  -  v)^,  B5  =  v‘ 
Bg  =  -2v(l  -  v),  =  5  -  6v,  C2  =  -(1  -  2v),  C3  =  2(1  -  v)  . 


The  strain-displacement  relations  as  shown  by  [27]  are 


2ej  •  =  w-  •  +  w.  .  g  W-  .  w,  . 

ij  1  .J  j,i  1,1  k,j 

with  ?  =  r  -  rQ  ,  ?  =  w  b^.  ,  w  =  g  ''w^ 


(2.8) 


where  b^  are  the  base  vectors  in  the  deformed  system,  and  r  and  r^  are 
the  position  vectors  in  the  deformed  and  undeformed  systems,  respec¬ 


tively. 


The  boundary  conditions  as  applied  to  the  stress-function,  $,  on 
a  free  surface  (e.g.,  along  the  crack  surface)  are  given  by  [28]: 


0  ,  $,2~0. 


(2.9) 


Thus  the  entire  system  of  equations  has  been  specified  for  an 
arbitrary  stress-type  problem  in  terms  of  orthogonal  curvilinear 


n 


coordinates  with  all  derivatives  taken  in  the  deformed  geometry. 

B.  Perturbation  Analysis 


The  problem  to  be  considered  is  shown  in  Figure  2.2,  and  consists 
of  a  small  ellipse  in  an  infinite  plate  with  all-around  uniform 
tension.  At  some  distance  from  the  ellipse  tip  the  stresses  are  re¬ 
duced  to  0(a),  the  strains  to  0(^/E),  and  the  displacements  to 
0(ca/E),  where  o'  is  the  remote  tension  and  c  is  one-half  the  crack 
length.  Therefore,  the  governing  equation  (2.7)  and  the  strain-dis¬ 
placement  relations  (2.8)  both  reduce  to  the  linear  forms.  Thus,  the 
solution  of  Inglis  [29]  applies 


,  sinh  25  (cosh  25  -  cosh  aSp) 
(cosh  2c  -  cos  2ri) 


(2.10a) 


sinh  25  .  (2.10b) 

(cosh  2c  -  cos  2n) 

.  .  ^  (cosh  2C  -  cosh  2Co) 
sinh  2n  TcosFtc  -  cbV"^  ’ 

{(1  -  2v)  cosh  2C  -  2(1  -  v)  cos  2n  +  cosh  2Co^  to 
e  - (cosh  2c  ■- 


(2.10c) 


w^  =  0  ,  (2.10e) 

where  C  and  n  are  the  elliptical  coordinates  defined  by  x  =  (cosh  c 
cos  n)»  y  =  sinh  c  sin  n;  e  =  ^(1  +  v),  the  characteristic  length  is  c, 
the  characteristic  stress  is  a,  and  Cq  is  the  non-dimensional  root 
radius  of  the  ellipsfe. 
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The  root  radius  of  a  crack  of  length  2c  after  deformation  is 

r  V2 

predicted  by  (2.10d)  with  Kq  =  0:w^  =  p  =  2e(l  -  v).  The  size  of 
the  ellipses  to  be  investigated  here  will  be  of  the  same  order, 
i.e.  5o  ~  0[2e(l  -  v)].  The  local  geometry  surrounding  the  tip  will 
then  be  as  shown  in  Figure  2.3  and  the  coordinates  will  locally  be 
reduced  to  parabolic: 


The  ordering  scheme  will  also  change  in  the  neighborhood  of  the  tip: 
the  stresses  are  0(1),  the  strains  are  0(1),  and  the  displacements  are 
0(1).  Therefore,  it  is  necessary  to  scale  the  variables  so  that  they 
are  all  0(1)  in  this  region: 


,  e^-j  ,  <\>  =  , 

“i  =  W,  ,  =  (2.11) 

Substituting  (2.11)  into  (2.7)  yields 


0  =  {1 

+ 

2Ai 

’nn 

2 

'^’’Cn 

+ 

A2 

(<}>>^^  +  <l>» 

^‘^’’CCCC 

+  2(f) 

’CCnn  ^  '^’’nnnri 

+ 

{1 

+  2(A3<|),^^ 

*  ^♦’nn 

)} 

2+ 

ccc 

®2*'Cnn' 

+ 

^3'*’ 

’CSC  “^’Cnn 

*  64*’nnn^^  %♦’ 

2  + 
CCn 

®6*’EEn*’nnn’ 

+ 

{1 

+  2(Ao((), 

'  3^’nn 

)} 

2+ 

nnn 

62*'SEn^ 

+ 

B3<}) 

’CCn'^’nnn 

%**EEE*’Enn’ 
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"  '=3(*-5nn*-nnn  *  ♦■5«*-E{„>  ’  (2''2) 


(The  notation  has  been  dropped  for  simplicity.) 
and,  the  substitution  of  (2.11)  into  (2.8)  leaves  it  virtually  un¬ 
changed. 

Since  the  displacements  in  the  tip  region  are  large,  the  boundary 
conditions  (i.e.,  ((>,•]  =  (j>,2  =  0)  must  be  applied  to  the  deformed 
geometry.  As  an  additional  boundary  condition  the  inner  solution  must 
be  'matched'  to  the  outer  solution  where  they  intersect.  As  Van  Dyke 
[30]  expresses,  the  outer  solution  (i.e.,  linear  solution  in  this  case) 
should  be  written  in  terms  of  the  inner  variables  and  then  the 
perturbation  parameter,  e,  should  be  made  very  small,  the  inner  solu¬ 
tion  should  then  be  matched  to  the  result.  The  inner  solution  would 
then  approach,  as  5  and  n  get  large,  these  expressions  for  stress, 
strain  and  displacement: 


;  ^  5(5^  +  2n2  + 

“  (C2  +  n2) 


(2.13a) 


.  _  r)U^  -  5o^) 

(?2+  r,2) 

5(5^  -  Co^) 

”5?  ■  (^2  +  ^2) 

(1  -  2v)  5^  +  2(1  -  v) 


(2.13b) 


{2.13c) 


(2.13d) 
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=  0  . 


(2.13e) 


As  the  problem  is  a  quasi -elliptic  fourth-order  partial  dif¬ 
ferential  equation,  there  must  be  two  boundary  conditions  at  each  of 
the  four  edges,  as  shown  in  Figure  2.4.  At  the  inner  edge  of  the 
ellipse  (5  =  Sg)  the  boundary  conditions  are  as  discussed  above: 

(|),g  =  0  and  (fi,^  =  0  or  alternately  <J),^  =  0  and  ^  =  constant  (where 
indicates  the  covariant  derivative  of  <J)  in  the  direction  normal 
to  the  ellipse  surface).  Since  the  problem  is  symmetric  about  the 
x-axis  the  boundary  conditions  along  n  =  0  reduces  to  |^  =  0  and 

=  0.  Now  the  matching  conditions  along  n  =  Ti,nax  and  5  =  must 
be  reduced  to  quasi -boundary  conditions.  Stress  boundary  conditions 
at  an  interface  between  two  materials  require  that  the  normal  normal 
stress  and  the  tangential  shear  stress  be  the  same  on  both  sides. 
Therefore,  these  conditions  shall  be  adopted  at  the  boundary  between 
the  inner  and  outer  regions  (referring  to  Figure  2.3): 

A  =  d,  d±  ^ 

linear  BS  ^^linear  ’  at  5  - 


M.  =  li 


linear  ’  "  9n,. 


at  n  =  T)r 


linear  , 


However,  for  matching  to  be  complete  the  inner  tangential  normal  stress 
at  the  interface  must  also  be  equal  to  the  outer  tangential  normal 


stress.  To  insure  this,  the  values  of  and 

iTIaX 

until  the  matching  occurs. 


will  be  enlarged 


Thus  the  solution  procedure  can  be  summarized  as  follows;  first, 
a  deformed  geometry  is  chosen;  then,  the  linear  outer  solution  is 
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computed;  next  the  non-linear  solution  is  determined  subject  to  the 
deformed  geometry  by  finite  differences  and  is  matched  to  the  outer 
solution  at  a  sufficiently  large  distance  from  the  crack  tip;  finally, 
the  non-linear  strain-displacement  relations  are  integrated  to 
determine  the  displacements. 


III.  SOLUTION  TO  THE  PROBLEM 

The  governing  equation  for  the  stress  function  will  now  be 
solved  for  various  ellipses  with  small  root  radii  approaching  that  of 
a  deformed  crack  in  an  infinite  plate.  In  order  for  the  nonlinear 
quasi -elliptic  equation  to  be  solved  it  will  be  first  reduced  to 
linear  form  by  a  Newtonian  method;  the  coefficients  of  this  linearized 
equation,  however,  will  depend  upon  the  stress  function  itself.  There¬ 
fore,  the  equation  must  be  solved  repeatedly,  updating  these  coef¬ 
ficients  each  time,  until  the  solution  changes  within  some  acceptably 
small  amount.  Since  the  convergence  of  the  equation  is  questionable 
and  because  of  the  lack  of  experimental  evidence,  two  general  solution 
methods  and  one  limit  analysis  will  be  presented.  The  first  method  to 
be  discussed  is  a  finite-difference  analysis;  the  second  is  an  inter¬ 
nal  collocation  analysis  based  upon  a  truncated  infinite  orthogonal 
series  which  satisfies  the  boundary  conditions  exactly;  and  finally, 
the  limit  analysis  which  is  based  upon  the  integration  of  the  strain- 
displacement  relations. 

A.  Linearization  of  the  Governing  Equation 

The  one  equation  for  the  stress  function  (2.12)  may  be 
represented  as 

0  =  AV^(j)  +  Ai  F]  +  A2F2  +  (3.1) 
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where 


A  =  1  +  2Ai(<}i,^^  +  (j),^^)  +  - 


+  A2(<|..55  +  ] 


A-,  =  1  +  2(A3<j),^^  +  ;  A2  =  1  +  2(A3<f,^^  +  A4<|),^^) 

Bi  2  Bp  2  ^3 


+  B4/B0  +  B5/B0  +  Bg/Bo  4>,g5^4>.^^^; 

F2  =  Bi/Bo  4^^^^  +  Bp/Bo  +  B3/B0 

+  B4/B0  +  B5/B0  +  Bg/Bp  4),55^4>,5^^; 


F3  =  Ci/Bo4,5^^4,^5^  +  C2/Bo4.5554>,^^^ 

+  ^3/Bo(4>.5Tin‘’’’nnn  ' 

Expanding  A  in  a  Taylor  series  about  some  initial  value  of  4>^^,  4>^^5 
and  4^^  yields 

A  =  Aq  +  ^^Aq(4,^^  -  4j^^q)  +  ^^^o('!’»5n  " 

*  -  ♦>nJo>  * 


where  ^  represents  the  initial  value 

and  6  =  (4,^5  -  4,5^13)  or  (4,^^  -  4,5^0)  or  (4,^^  -  4,Jo) 

with  =  14;^=  2A^  +  A(4,^^  +  2A2(4,g5  +  4,^^)) 

- =  -84,. _ 
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In  like  manner  linearizing  substituting  these  into 

(3.1)  yields 


+  ^6  (b,^^  +  °3  (|),r  +  '6  (t,  =  6 

0  55  o^’5n  0  in 

where  o$  =  A;  h  = 


(3.2) 


^3  =  +  ??^F2A2  +  ^^'^F3A3;  h  =  ^"^^F^Ai  +  ^'^’^F2A2 

+  ^^'^F3A3;  ^3  =  ^"^lA.,  +  ^^’1F2A2  +  '^’''^F3A3; 

^3  =  ^^A.|Fi  +  ^^A2F2  +  ^^AV^(f.;  ^3  =  ^’^A3F3  +  ^'^AV^<() ; 

^3  =  ’^’^AiF]  +  '^'^A2F2  +  '^'^AV^({);  ^3  “  ”(^1^1  ^  ^2^2  ^  ^3^3^o 
*  'V’555o  ^  ^®o*-ECno  *  ^'’o^’Enno  ^  '’^o+’n™ 

^  V-Eno  ■"  V-nno’ 


(A  complete  listing  of  the  quantities  etc.,  and  the  necessary 

physical  covariant  derivatives  are  listed  in  the  Appendix.) 

B.  Finite-Difference  Method 


The  governing  equation,  having  now  been  reduced  to  a  tractable 
linear  form,  can  now  be  solved  by  finite  differences.  Near  the  boundary 
of  the  nonlinear  region  (the  boundary  comprising  the  inner  stress-free 
border  at  or  the  linear-nonlinear  interface)  the  finite  difference 
scheme,  written  at  point  2  as  shown  in  Figure  3.1a,  becomes 
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m  ='^i^+o(Ax2) 

dx  2AX  '  ^ 


Various  two-dimensional  combinations  of  the  two  above  mentioned 
schemes  are  used  depending  upon  the  placement  of  the  grid  point  (at 
which  the  governing  equation  is  being  solved)  in  relation  to  the 
boundaries.  Two  specific  cases  of  those  possible  will  be  discussed 
here:  first,  the  scheme  used  when  no  boundaries  are  near;  and  second, 
the  scheme  used  when  two  boundaries  are  near.  In  the  first  case  as 
shown  in  Figure  3.2a  the  equation  will  be  solved  at  point  13  with  the 
resulting  expressions  for  the  partial  derivatives: 


3^1  _  (<t’15  ~  ^'{’14  ^*^13  ~  ^^12  ‘i^ll)  +  o(.2) 

Av4 


4^ _  _  (<i>i9  -  2<()-|3  +  ^■\y  -  2((.-]4  +  4(|)-|3  -  2<i)]2  +  <l>9  - 

2  2  A  A  \/2 


3x^3y 


+  0(62) 


8  ^^13 

A  \/4 


^  0(62) 


(<)>15  -  2(j)i4  +  2(^■^2  -  <i>i] ) 


+  0(6^) 


3x^3y 


^  (<^19  -  2j.^8  +  <I>17  -^q  +  24q  -  ^7)  ^  ^^^2^ 


2Ax^Ay 
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9xay'^J^3 


19  ~  ‘1*17,  ~  ^*^14  ^‘*’12  '*'  *^9  ~  <t'7^  ^  (3  5) 


2AxAy 


a^d) 


9y'' 


=  (<^23  ~  ^^18  +  ^‘i>8  ~  *1*3)^^  0(^2) 

13 


ax^j 


13 


=  (<^14  ~  ^*<>13  ^  '^12)  +  0(6^) 

AX^ 


a^d)  =  (<j>l9  -  <t>l7  -  <^9  +  ‘t’j)  ^  q/52x 
axay  4AxAy  '  ^ 


13 


a^d)'' 


ay^J 


_  (<^18  -  ^<<>13  +  ^s)  ^  q(^2) 
13 


M. 

ax 


(<f>14  -  '}>-|2^ 


2ax 


+  0(6^) 


13 


ii 

ay 


13 


,  ,<*18  :.*$l+0(s2) 

2Ay 


vjhere  S  is  either  ax  or  Ay. 


For  the  second  case  when  two  boundaries  are  present  as  shown  in  Figure 
3.2b,  the  governing  equation  will  be  solved  at  point  9  with  the  follow 
ing  expressions  for  the  partial  derivatives  resulting: 


a'^d> 


(-113d>iQ  +  60Ax||-)^q  +  192d)9  ~  108d>3  +  324>y  - 

12ax^ 


+  0(6^) 
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3^(J) 


(4'15  "  2(|)i^  +  (^12  -  2i})iq  +  44)g  -  2<t)Q  +  '{>5  -  2<^4  +  ({>3) 


2  2 
&x  Ay^ 


0(62) 


a^cf. 


sy^J 


(-34'24  +  32^19  -  108(|)]4  +  192(j)9  -  60Ay|^)^  -  n3<l>4)  ^  2^ 

12Ay4  ^  ^ 


dx- 


*  *  ^*9  -  *7)  ^  „(j2) 

3ax2 


ax2ayj 


3^4) 


(4>15  -  24)14  +  <i>13  '  +  2(|)4  -  4)3)  ^  ^^^2) 


2Ax^Ay 


3x3y2j 


=  («l’15  ~  *<’13  '  2<^io  +  24>g  +  4>5  -  <(>3)  ^  q^^2) 
2AxAy2 


(3.6) 


3^4» 


3y' 


(<}»19  -  94>9  +  6Ay|^)^  +  84.4  ^ 
3Ay2 


3^4) 


3X‘ 


iMJL^^aJLial  +  o(«2) 

AX^ 


324) 


3x3y 


(4,15  -  4>13  -  <l>5  +  <l>3)  ^  2x 

4AxAy  ' 


3^]  ^  (*14  -  ^*9  *  *4)  ^  ^jj2j 


ay^ 


Ay‘ 
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J  =  t  0(s2) 

3X  2AX 

^9 


J  ,  .(.♦.14  -  ■>4)  ^  o(;2)  . 

9y  2Ay 

Jg 

\ 

M. 

(Note  that  the  values  of  (^>2^9  4>20^  ^15»  ^10’  ^5’  ^4’  ^3^  ^2*  ’  9x 

>10 

\ 

and  are  all  known  and  treated  as  boundary  conditions.) 
dyj4 

Using  equations  (3.5)  or  (3.6)  or  an  anologue  of  them,  all  of 
the  partial  derivatives  necessary  in  the  evaluation  of  the  needed 
physical  covariant  derivatives  in  the  governing  equation  (3.2)  can  be 
obtained.  The  region  of  interest  should  be  divided  up  evenly  in  the 
parabolic  coordinates  ^  and  n  to  give  the  necessary  grid  pattern  as 
shown  in  Figure  3.3.  Then  at  each  grid  point  interior  to  the  boundary 
a  linear  equation  involving  a  bandwidth  of  3N  +  1  (where  N  is  the 
number  of  grid  points  in  the  n  direction)  is  generated.  This  resulting 
linear  system  of  equations  may  be  solved  simultaneously  by  any  direct 
method,  e.g.  Gaussian  elimination: 

CA](t  =  B  (3.7) 

where  [A]  is  a  banded  coefficient  matrix  depending  upon  the  distri¬ 
bution;  (f)  is  the  independent  variable  evaluated  at  each  grid  point; 

B  is  a  vector  of  constants  composed  of  in  equation  (3.2)  and  the 
boundary  conditions.  Since  the  coefficient  matrix  and  the  right-hand- 
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side  vector  both  depend  upon  the  ij)  distribution,  the  solution  must  be 
iterated  until  the  change  is  acceptably  small. 


C.  Internal  Collocation  Method 


The  second  method  pursued  involves  selecting  a  series  of 
orthogonal  polynomials  which  satisfy  the  boundary  conditions  exactly. 
One  such  representation  is  the  familiar  double  Fourier  series: 


^  ~  “^linear 


+ 


M 

I 

m=l 


N 

E 

n=l 


"mn 


{1  -  cos 


{1 

'^max  ”  ^0^ 


(3.8) 


Note  that  this  double  series  satisfies  the  symmetry  condition  at  n  =  0, 
the  boundary  conditions  at  ?  =  Cq*  the  matching  conditions  at 
K  =  5max>  the  matching  conditions  at  n  =  nmax-  There  are  MxN 
unknown  constants  in  this  double  series;  these  may  be  determined  by 
solving  the  governing  equation  (3.2)  at  MxN  discrete  points  in  the 
region  of  interest  by  iteration.  Figure  3.3  shows  how  the  region  is 
to  be  divided  up  evenly  (equal  divisions  are  not  necessary;  hov/ever, 
convergence  of  the  method  were  helped  by  this).  The  resulting  system 
of  linear  equations  may  be  solved  by  direct  methods: 

[A]Cfn„  =  B  (3.9) 

where  A  is  a  full  MxN  by  MxN  coefficient  matrix  depending  upon  the 
distribution;  C^p  are  the  coefficients  in  equation  (3.8);  and  B  is  the 
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vector  of  constants  depending  on  in  the  previous  iteration. 

The  advantage  of  this  method  lies  in  its  simplicity  of  formu¬ 
lation  and  in  its  satisfaction  of  the  boundary  conditions  exactly. 
However,  due  to  the  ripples  inherent  in  any  Fourier  representation, 
the  method  became  unstable  for  large  values  of  M  and  N.  Also,  since 
the  coefficient  matrix  [A]  in  equation  (3.9)  is  completely  full,  the 
solution  time  and  core  size  became  large  for  even  moderate  numbers  of 
terms.  The  finite-difference  formulation  is  complex  because  of  the 
many  types  of  grids  possible;  however,  the  method  is  stable  and  the 
coefficient  matrix  [A]  in  equation  (3.7)  is  sparse  and  can  be  solved  by 
banded  Gaussian  methods  with  a  considerable  savings  in  computing  time 
on  a  digital  computer. 

D.  Integration  of  the  Strain-Displacement  Relations 
and  Limit  Analysis 


Whichever  method  is  used  to  determine  the  stress  function  and 
subsequently  the  covariant  strains  e^j,  these  strains  will  be  inte¬ 
grated  through  the  complete  strain -displacement  equations  (2.8)  to 

I 

determine  the  displacements.  Writing  equation  (2.8)  in  terms  of  the 
parabolic  coordinates  yields 


2  2 

w^,^/g  - 


(3.9a) 


Hr,  “  «{-n  *  “n’S  -  -  "n’n'''n-E^« 


(3.9b) 

(3.9c) 


(The  notation  has  been  dropped  for  simplicity.) 
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These  equations  (3.9)  will  be  finite-differenced  to  yield  a  system  of 
first  order  nonlinear  differential  equations  which  must  be  solved 
simultaneously  through  iteration.  The  finite  difference  scheme  to  be 
used  here  is  shown  in  Figure  3.4  and  the  resulting  partial  derivatives 
solved  at  point  3  are 


Since  there  are  only  two  unknowns  (w^  and  w^)  to  be  determined,  only 

two  of  the  three  equations  (3.9)  are  necessary  in  the  procedure. 

Equations  (3.9a)  and  (3.9c)  are  chosen  because  of  their  superior 

stability  in  the  iterative  calculations.  Solving  (3.9a)  for  w 

and  (3.9c)  for  w  yields 
n’n 

■  9  (1  *  •'l  -  w^25/g2  -  2£;;/g  ) 


w 


g  (1  ± 


The  lower  sign  will  be  chosen  because  the  expression  would  then  reduce 
to  the  linear  result  for  small  strains  and  displacements;  hence,  the 
above  result  would  become 
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w„,„/g  =  (1  -  •'1  -  ) 


Introduction  of  a  to  represent  a  physical  quantity  gives  the  follow¬ 
ing  expressions: 


“c.f  '  <'  ■  > 


V/  =  (1  -  *^1  -  W  ^  -  2e  ) 

n,n  '  ^>r\  nn 


(3.11a) 

(3.11b) 


(The  necessary  physical  covariant  derivatives  of  w^  and  w  are  found 

C  n 

in  the  Appendix.) 

After  substituting  the  finite  difference  equations  (3.10)  into 
the  linearized  expressions  for  and  in  equations  (3.11),  the 
actual  numerical  procedure  to  find  the  displacements  would  then  begin. 
Referring  to  Figure  3.3,  the  integration  procedure  v/ould  begin  one 
step  from  Cmax  ~  5max  ~  at  n  =  0.  Since  at  this  grid  point 

w^  =  0  equals  zero  and  w^^^  =  0  equals  zero  from  symmetry  conditions, 
equation  (3.11a)  becomes  uncoupled  and  can  be  solved  for  w^  directly 
without  iteration,  j Next  the  procedure  steps  in  the  positive  rrdirection 
solving  equations  (3.11a)  and  (3.11b)  simultaneously  for  w^  and  w^ 
through  iteration  at  each  grid  point  until  rimax  reached.  Then 
another  step  is  taken  toward  Cq  (Snew  ""  Sold  "  ^s)  and  the  procedure 
is  repeated.  This  stepping  in  both  the  negative  g  and  the  positive 
n-di recti ons  is  carried  out  until  w^  and  w^  have  been  determined  in 
the  entire  region. 

The  limit  analysis  consists  of  examining  the  existence  of 
in  equation  (3.11b)  at  the  ellipse  tip  along  the  line  of  symnetry. 
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~  2  ~ 

The  discriminant  of  this  equation  is  1  -  Noting  that 

is  equal  to  zero  along  the  line  of  symmetry  reduces  the  expression 
to  1  -  2e^^.  From  equation  (2.4)  can  be  written  in  terms  of  the 
stresses: 

^nn  =  -  ^)  %n  ' 

Along  the  ellipse  equals  zero,  so  that  -  (1  -  v)  0^^. 

Inserting  this  result  into  the  expression  for  the  discriminant  above 
yields 

discriminant  =  1  -  2(1  -  v)  0^^. 

For  to  be  real  in  equation  (3.11b)  the  discriminant  must  be 
greater  than  or  equal  to  zero: 


0^1  “2  (1  “v)  9t  ^  ^0’  ~  ^ ’ 

or  alternately  (3.12) 

0  <  oTT^ - T  3t  5  =  5g,  n  =  0. 

Tin  -  2(1  -  v)  ^ 


To  understand  the  result  expressed  in  equation  (3.12)  more  fully 
it  is  necessary  to  study  the  unstrained  metric  tensor  h^j  from 
equation  (2.2b): 


hlj  =  9 


0  -  2e„) 


-2e 


Sn 


The  determinant  of  this  matrix  may  be  written 
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|h|  =  g2{(l  -  2?jj)  (1  -  2e^„)  -  455/}; 

along  the  line  of  symmetry  this  expression  may  be  simplified  to 

|hl  =  g2{{l  -  2e^^)  (1  -  2e^^)}  (3.13) 

If  the  determinant  |h|  is  zero,  this  situation  is  equivalent  to  a 
mapping  which  takes  a  region  in  the  deformed  geometry  and  transforms 
it  into  a  point  in  the  initial  geometry,  i.e.  an  ellipse  maps  into  a 
crack  (a  degenerate  ellipse).  At  the  ellipse  tip  becomes  zero  so 
that  equation  (2.4)  would  give  the  strains: 

"nn  =  -  ^^%n  ’ 

substituting  this  result  into  (3.13)  and  equating  this  result  to  zero 
gives 

0  =  (1  +  v5^^)  (1  -  2(1  -  v)5^^)  . 

Because  the  plate  is  under  uniform  tension  it  is  unlikely  that 
would  ever  be  negative,  the  above  equation  becomes 

0  =  (1  -  2(1  -  v)5^^)  or  •  (3.14) 

Comparing  the  results  of  (3.14)  and  (3.12)  it  is  evident  that  a  value 
of  higher  than  is  impossible  because  the  ellipse  would  map 

past  a  crack.  Therefore  the  maximum  stress  possible  at  the  tip  of  an 
ellipse  is  "  2(1  -  v')’  occurs  if  the  initial  geometry 

is  a  crack  of  zero  root  radius. 


IV.  DISCUSSION  OF  THE  RESULTS  AND  CONCLUSIONS 

The  physical  stresses  and  strains  in  the  field  surrounding  the 
tip  will  now  be  presented  for  small  ellipses  with  Poisson's  ratio 
equal  to  one-half;  the  contravariant  displacements  in  the  tip  region 
will  also  be  discussed.  Finally,  the  results  for  these  small  ellipses 
will  be  extrapolated  to  determine  the  stress  field  around  the  deformed 
tip  of  a  crack  of  zero  initial  root  radius. 

Throughout  the  following  discussion  reference  will  be  made  to 

the  terms  'physical'  and  'dimensional'.  A  'physical  dimensional' 

stress  will  be  a  stress  as  measured  experimentally  and  will  be  denoted 

a.{j.  A  'physical  nondimensional '  stress  can  be  found:  aij  =  o^-j/a, 

where  is  the  remote  'physical  dimensional'  stress.  A  covariant 

stress  CTij  may  be  written  as  +  n^.  Finally  the  scaled 

inner  variables,  e.g.  ^ij,  will  not  have  the  star  superscripts  (*) 

★ 

whereas  the  outer  variables  will  have  them,  e.g.  o-jj: 

r  w5*  pj?l 

^ij  =  ’  "ij  =  ^ij*  ’  w  =  e  ’  P  -  e2  ’ 

y*  y*  o 

X  =  -r  ;  y  =  2  >  ^  r 

e 

A.  Discussion  of  the  Results 

The  selected  root  radii  of  the  ellipses  in  the  deformed  state 
range  from  p  =  1  to  9  which  corresponds  to  a  Sq  1  to  3.  The  minimum 
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root  radius  (p  =  1)  is  that  of  a  deformed  crack  as  predicted  by 
Inglis  [29].  The  physical  normal  stress  distribution  along  the 
line  of  symmetry  in  front  of  the  tip  for  =  1.0,  1.3,  1.7,  and  2.5 
is  shown  in  Figures  4.1a,  4.2a,  4.3a,  4.4a.*  The  linear  curve  referred 
to  here  and  afterwards  concerns  the  analysis  of  Inglis.  The  physical 
normal  stress  distribution  along  the  ellipse  border  for  Cq  “ 

1.3,  1.7  and  2.5  is  shown  in  Figures  4.1b,  4.2b,  4.3b  and  4.4b. 
Examination  of  Figures  4.1a  and  4.1b  show  the  relative  merit  of  the 
finite  difference  as  opposed  to  the  collocation  methods  for  solving 
the  nonlinear  governing  equation.  As  mentioned  previously  in  Chapter 
III,  for  a  small  number  of  terms  the  collocation  approach  gave  fair 
results;  however,  for  larger  numbers  of  terms  the  approach  became  un¬ 
stable.  Therefore  the  finite  difference  approach  is  viewed  as  being 
more  accurate  although  they  both  give  qualitatively  the  same  results. 
Referring  to  Figure  4.1a  with  linear  solution  falls 

monotonically  from  2.0  at  the  tip  of  the  deformed  ellipse.  The  non¬ 
linear  (finite  difference)  curve  rises  from  very  close  to  1.0  to  a 
maximum  of  1.24  and^then  falls,  approaching  the  linear  curve  from  the 
top.  The  maximum  error  then  is  50%  and  this  occurs  at  The  limit 
analysis  predicts  that  the  value  of  at  the  ellipse  tip  should  be 
less  than  or  equal  to  1.0  and  this  agrees  with  the  nonlinear  results 
to  within  1%.  Referring  to  Figure  4.1b  with  ^ linear 
and  the  nonlinear  curves  fall  monotonically  with  the  linear  curve  being 
always  above. 

*In  all  of  Figures  4,  points  where  inner  and  outer  solutions 
were  matched  fell  well  outside  regions  shown. 
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The  physical  normal  stress  distribution  along  the  line  of 
symmetry  in  front  of  the  tip  for  =  1.0,  1.3,  1.7  and  2.5  is  shown 
in  Figures  4.1c,  4.2c,  4.3c  and  4.4c.  The  disparity  between  the  linear 
and  the  nonlinear  analysis  for  the  physical  normal  stress  with 
=  1.0  along  the  line  of  symmetry  can  be  seen  in  Figure  4.1c.  The 
maximum  difference  between  the  two  distributions  is  0.05  and  is  not 
significant;  although,  the  disparity  represents  a  large  percentage 
because  of  the  small  magnitude  of  the  stress. 

In  linear  fracture  mechanics  the  single  parameter  of  interest 
is  the  stress  intensity  factor  K  which  defines  the  field  strength  of 
the  singular  behavior.  For  a  crack,  as  shown  in  Figure  2.2,  K  as 
defined  by  Sneddon  and  Lowengrub  [31]  is  lim*^2x*  a  (x,0).  The  linear 

X-K)_ 

solution  with  5o  =  0  would  then  give  as  •  Defining  an 
"apparent"  stress  intensity  factor  as  Ky^p  =  *'2>r*  ^  »  then  one  obtains 
Kyvp 

■j^  =  /2>r  Hence,  in  Figures  4. Id,  4. 2d,  4.3d  and  4.4d  for 

Cq  =  1.0,  1.3,  1.7  and  2.5  the  apparent  stress  intensity  factor  K/\p 
is  shown,  with  the  'singular'  curve  representing  the  linear  solution' 
to  the  undeformed  crack  geometry  (Sq  =  0)-  The  error  involved  in  the 
calculation  of  stress  intensity  factors  for  Sq  “  ^  ''S  shown  in  Figure 
4.1d.  The  nonlinear  curve  represents  the  K/\p  that  would  be  calculated 
from  data  taken  from  an  actual  experiment.  If  the  data  were  taken  out¬ 
side  of  X  =  16,  there  would  be  only  an  error  of  2%.  However,  the  error 
would  rise  to  +20%  and  then  drop  to  -100%  as  x  decreases  to  zero.  The 
linear  analysis  predicts  that  the  error  would  rise  to  +10%  before 
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dropping  to  -100%.  These  results  are  consistent  v/ith  recent  experi¬ 
mental  findings  [10],  [14],  [18]. 

The  physical  maximum  in-plane  shear  stress  (the  primary  quantity 
used  in  photoelastic  experiments)  along  a  line  originating  at 
the  tip,  perpendicular  to  the  line  of  symmetry,  for  Cq  =  1-3, 

1.7  and  2.5  is  shown  in  Figures  4.1e,  4.2e,  4.3e  and  4.4e.  The  dif¬ 
ference  between  the  nonlinear  and  the  linear  formulations  with  ?o  “  ^ 
for  is  shorn  in  Figure  4.1e.  The  maximum  difference  between  the 
linear  and  the  nonlinear  is  +50%  and  occurs  at  the  tip.  Therefore, 
the  linear  approach  represents  the  distribution  v;ell  until 
y  <  l~p;  i.e.  inside  of  one  root  radius  it  is  necessary  to  go  to  the 
nonlinear  representation. 

After  integrating  the  strain  displacement  equations  (3.9)  for 
the  covariant  displacements  (w^,w^)  the  contravariant  counterparts  can 
then  be  determined: 


w^ 


wn 


These  displacements  can  then  be  used  to  find  the  initial  geometry: 

(4.1a) 

-  w*^  ,  (4.1b) 


^initial  “  ^final  " 


^initial  ~  ^final 


The  contravariant  displacement  w^  distribution  along  the  line  of  sym¬ 
metry  for  5o  =  1-3,  1.7  and  2.5  is  shown  in  Figures  4. If,  4.2f, 

4.3f  and  4.4f.  The  difference  between  the  linear  and  the  nonlinear 
curves  for  Sq  “  ^  shown  in  Figure  4. If.  The  error  is  within  3% 
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until  very  close  to  the  tip:  here  the  error  rises  to  8%.  The  linear 
displacements  would  map  this  ellipse  exactly  into  a  crack;  and  there¬ 
fore,  the  small  difference  between  the  displacement  curves  indicates 
that  the  nonlinear  analysis  of  an  ellipse  of  Sq  ”  ^  would  closely 
resemble  that  of  a  crack  (in  the  initial  geometry),  except  near  the 

tip  itself. 

Figure  4.5  shows  the  variation  of  the  physical  normal  stress 
distribution  along  the  line  of  symmetry  for  the  range  of  final, 
deformed  root  radii  p  =  1-^9.  From  this  figure  it  can  be  seen  that 
the  maximum  value  of  the  physical  normal  stress  occurs  at  the 
ellipse  tip  (5  =  Cq)  for  5o  >  ‘■•'7’  however,  for  Co  <  the  maximum 
occurs  at  some  distance  (C  >  Cq)  phenomenon 

suggests  two  different  types  of  crack  growth:  1)  if  the  root  radius 
of  the  notch  is  above  some  critical  level,  the  flaw  will  extend  from 
the  edge  of  the  blunted  notch  as  reported  by  Schijve  and  Jacobs  [32]; 
2)  if  the  root  radius  of  the  notch  is  below  some  critical  level,  the 
flaw  will  extend  by  the  creation  of  crack  segments  ahead  of  the  sharp 

notch  as  cited  by  Yokobori  [33]. 

The  distribution  of  the  physical  normal  strain  along  the 
line  of  symmetry  in  front  of  the  tip  for  the  range  of  final  root  radii 
p  =  1^9  is  shown  in  Figure  4.6.  The  limit  analysis  showed  that  a 
physical  strain  value  of  0.5  was  impossible;  in  fact,  this  value 
would  only  be  reached  near  the  crack  tip.  The  above  figure  bears  out 
this  conclusion  in  that  all  of  the  curves  are  asymptotic  to  0.5  at 

K  -  ^0’ 
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From  the  contravariant  displacements  shown  in  Figures  4.1 f 
through  4.4f  the  initial  geometry  of  the  ellipses  can  be  determined 
through  equations  (4.1).  Via  this  procedure  the  initial  root  radii 
of  the  deformed  ellipses  are  found  and  they  are  presented  in  Figure  4.7. 
The  linear  analysis  as  presented  in  this  figure  would  give  an  ellipse 
as  the  initial  geometry;  however,  the  nonlinear  analysis  would  not 
and  the  difference  between  the  two  curves  represent  the  amount  the 
initial  geometry  would  differ  from  an  ellipse.  From  the  extrapolation 
of  this  curve  it  is  clear  that  a  crack  initially  with  no  load  would 
deform  into  a  quasi-elliptical  shape  with  a  root  radius  of  (0.92)^. 

From  the  displacements  for  ?o  “  1*0  it  is  clear  that  the  crack  will 
deform  into  an  ellipse  of  shape  =  1.0  at  some  distance  from  the  tip. 
It  can  be  stated,  therefore,  that  the  behavior  in  the  neighborhood  of 
the  tip  can  be  well  approximated  by  an  ellipse  of  Kq  -  0.92  and  that 
away  from  the  tip  the  behavior  can  be  determined  by  an  ellipse  of 
Co  “  1  •  0 . 

Using  Figures  4.5  and  4.7  the  magnitude  and  the  location  of  the 

maximum  normal  stress  can  be  plotted  with  respect  to  the  initial 

^%ax 

root  radius  as  is  shown  in  Figures  4.8  and  4.9,  respectively.  Ex¬ 
trapolating  the  curve  on  Figure  4.8  the  maximum  normal  stress 
becomes  1.3,  approximately;  and  referring  to  Figure  4.9  the  location 
of  the  maximum  stress  becomes  x/p  =  0.33,  so  that  x  =  0.33(0.92)^  = 

0.28.  It  could  be  conjectured  that  the  shift  of  the  maximum  stress  in 
front  of  the  tip  could  be  as  a  result  of  the  additional  constraint 
in  the  interior  of  the  body. 
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The  near  field  zone  of  the  present  model  extends  to  O(e^)  from  the 
crack  tip  as  contrasted  with  0(e)  in  the  Rice-Rosengren  [21]  and  Hutch¬ 
inson  [22]  models,  and  therefore  is  a  very  near  field  solution.  The 
model  reveals  that  the  strains  are  limited  as  shown  in  Figure  4.6; 
physically  this  means  that,  upon  application  of  small  remote  stresses, 
the  portion  of  the  crack  surfaces  nearest  the  tip  each  undergo  a  rota¬ 
tion  of  90°  in  opposite  directions  in  creating  the  notch  tip.  This 
produces  stresses  of  the  order  of  Young's  modulus  at  the  notch  tip  (Fig.  4.5) 
as  in  [25]J26]  but  these  local  stresses  are  independent  of  the  applied 
remote  stress.  Increasing  the  remote  stress  simply  enlarges  the  root 
radius,  and  with  it  the  size  of  the  zone  affected  by  the  finite  defor¬ 
mations.  These  features  (i.e.  stresses  of  the  order  of  the  modulus 
and  independence  of  the  near  field  stress  from  the  remote  loads)  are 
also  found  in  the  "cohesive"  stress  field  of  the  Barenblatt  model  [34]. 
Moreover,  the  "hump"  in  the  a  stress  distribution  shown  in  Figure 
4.5for  small  root  radii  produces  qualitatively  the  same  sort  of 
stress  distribution  very  near  the  notch  tip  as  Neuber  obtained  with 
his  "limiting  particle  size"  theory  for  pointed  notches  [35]. 

For  stress  freezing  photoelastic  materials  above  critical  tem¬ 
perature,  the  maximum  stress  is  expected  to  be  one  order  of  magnitude 
below  Young's  modulus  and  this  adjustment  might  possibly  be  incor¬ 
porated  into  the  present  model  by  introducing  non-linear  constitutive 
relations.  This  change  would  also  be  expected  to  enlarge  the  size  of 
the  affected  zone  (currently  of  the  order  of  one  root  radius)  to 
perhaps  several  root  radii. 
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B.  Conclusions 

A  finite  deformation  analysis  of  the  near  field  surrounding  the 
tip  of  small  crack-like  ellipses  was  performed  in  the  deformed  geometry 
using  the  complete  compatibility  equations  and  linear  constitutive 
relations. 

The  effects  of  finite  strains  and  rotations  in  the  tip  region 
of  a  deformed  crack  were  found 

(1)  to  reduce  and  at  the  tip  by  50%; 

(2)  to  move  ^Tinp^ax  " 

1.3); 

(3)  to  reduce  pfinal  (0.92)  ; 

(4)  to  change  insignificantly. 

These  contributions  are  all  extremely  local  and  are  concentrated  within 
a  few  root  radii  of  the  tip  and  would  not  normally  be  measurable  in 
a  photoelastic  experiment;  however,  the  effects  of  crack  tip 
blunting  are  more  general  and  spread  back  to  ten  or  more  root  radii. 

For  initial  geometries  with  blunted  tips,  the  regions  of  influence 
cited  above  are  generally  the  same  (in  terms  of  root  radii)  although 
the  magnitudes  of  the  above  effects  are  less.  The  above  results 
suggest  that  blunt  notches  propagate  from  the  notch  tip;  however, 
sharp  notches  extend  from  the  formation  and  subsequent  growth  of 
small  cracks  in  front  of  the  notch. 

The  features  of  the  present  model  are  believed  to  be  qualitatively 
appropriate  to  behavior  of  photoelastic  materials  above  critical  tempera¬ 
ture.  Studies  are  currently  underway  for  the  inclusion  of  non-linear  con¬ 
stitutive  relations  in  order  to  improve  its  quantitative  predictions  and 
hopefully  shed  new  light  on  the  large  strain  region  very  near  the  crack  tip. 
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APPENDIX  A 


The  coefficients  as  they  appear  in  equation  (3.2)  are 
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for  the  elliptical  coordinates 
X  »  cosh  C  coshiy  ,  y=  sinhf  sinhi? 

°0*h  -'=°« 
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Figure  2.2  Problem  Geometry 


Figure  2.3  Inner  Region  Geometry 
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(a)  Scheme  Used  Near  a  Boundary 


(b)  Scheme  Used  Away  From  a  Boundary 

Figure  3.1  One-Dimensional  Finite  Difference  Schemes 
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(a)  Scheme  Used  Away  From  All  Boundaries 


Figure  3.2  Two-Dimensional  Finite  Difference  Schemes 
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Figure  3.3  Grid  Pattern 
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Figure  4.1a  Physical  Normal  Stress  a  Distribution  Along  the  Line  of  Symmetry  for  ^ 


Figure  4.1b  Physical  Normal  Stress  a  _  Distribution  Along  the  Ellipse  Border  for 


0.5 


Stress  Cpr  Distribution  Along  the  Line  of  Symmetry 


NON-LINEAR 


Figure  4.1e  Physical  Maximum  In-Plane  Shear  Stress  Distribution  for  Cq 


Stress  a  Distribution  Along  the  Line  of  Symmetry 


Figure  4.2b  Physical  Normal  Stress  a  Distribution  Along  the  Ellipse  Border  for 


Figure  4.2c  Physical  Normal  Stress  o  Distribution  Along  the  Line  of  Symmetry  for 


NON -LINEAR 


Apparent  Stress  Intensity  Factor  Distribution  Along  the  Line  of  Symmetry 


Figure  4.2e  Physical  Maximum  In-Pla 


NON-LINEAR 


Figure  4.2f  Contravariant  Displacement  w^Distributi 


T 


Figure  4.3a  Physical  Stress  a  Distribution  Along  the  Line  of  Symmetry  for 


Distribution  Along  the  Ellipse  Border  for 


Stress  o,.  Distribution  Along  the  Line  of  Symmetry  for 


NON-LINEAR 
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Figure  4.3d  Apparent  Stress  Intensity  Factor  Distribution  Along  the  Line  of  Symmetry  for  C 


LINEAR 


Contravan'ant  Displacement  w  Distribution  Along  the  Line  of  Symmetry 


Figure  4.4a  Physical  Normal  Stress  a  Distribution  Along  the  Line  of  Symmetry  for 


X) 


Figure  4.4b  Physical  Normal  Stress  5  Distribution  Along  the  Ellipse  Border  fo 


Stress  Oj-i-  Distribution  Along  the  Line  of  Symmetry 
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Figure  4.4e  Physical  Maximum  In-Plane  Shear  Stress  Distribution  for  ^0 


linear 


Figure  4.4f  Contravariant  Displacement  Distribution  Along  the  Line  of  Symmetry  for 


X//D  qt  (Tr^rj  =  (^7777)|viax 
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